Vibrationally Induced Two-Level Systems in Single-Molecule Junctions 
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Single-molecule junctions are found to show anomalous spikes in dl/dV spectra. The position in 
energy of the spikes are related to local vibration mode energies. A model of vibrationally induced 
two-level systems reproduces the data very well. This mechanism is expected to be quite general for 
single-molecule junctions. It acts as an intrinsic amplification mechanism for local vibration mode 
features and may be exploited as a new spectroscopic tool. 
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A single atom or molecule with an almost transparent 
single conductance channel leading to a conductance near 
the conductance quantum 2e 2 /h (= 1 Go) can be con- 
tacted to leads. Conduction electrons can pass through 
such junction ballistically for low bias voltages since the 
mean free path of the electrons is much larger than the 
size of the contact. However the contact is not entirely 
ballistic in the sense that once the excess energy of the 
conduction electrons becomes equal or larger than the en- 
ergy of a local mode of the contact, the electrons can scat- 
ter inelastically by exciting a local mode. This results in 
the case of a perfectly transmitting single channel contact 
to a small decrease in the conductance, since the forward 
travelling electrons are backscattered due to the energy 
loss in the inelastic scattering process. Differential con- 
ductance (dl/dV) measurements have identified vibra- 
tion modes of single molecules in an atomic contact Q. 
This technique, also called Point Contact Spectroscopy 
(PCS) is analogous to inelastic electron tunnelling spec- 
troscopy (IETS) for single molecules 0,0], with the dif- 
ference that the conductance in the latter case increases 
due to the opening of an additional conductance channel. 

In this letter we present the observation of anomalous 
spikes, rather than steps, in dl/dV measurements on var- 
ious single molecule contacts. We present a model that 
involves two-level systems, which describes our data very 
well. It may be used as a new spectroscopic tool for 
identifying molecular vibration modes in single molecule 
junctions. 

We create atomic contacts using a mechanically con- 
trolled break junction (MCBJ) setup in cryogenic vac- 
uum at 4.2 K (see ref. for a detailed description). 
Break junctions for the metals Au, Ag, Pt and Ni have 
been investigated with the molecules H 2 , D 2 , O2, C2H2, 
CO, H2O and benzene. In most of these cases regular vi- 
bration mode spectra displaying a step down in conduc- 
tance have been observed, but for all systems anomalous 
spectral features as displayed in Fig. were also found. 

In order to admit these molecules to the metal atomic 
contacts at 4.2 K, the insert is equipped with a capil- 
lary that has a heating wire running all along its interior 
to prevent premature condensation of the gasses. The 
amount of gas admitted is of order 10 /imol. Previous 



measurements have clearly demonstrated that it is pos- 
sible to capture a single molecule in the atomic junction 
and measure its vibration modes 0, . 

Figure QJa) shows a dl/dV spectrum of a single CO 
molecule contacted by Pt leads. The observation of spikes 
that dominate the spectrum is surprising. Pt-CO-Pt 
contacts at a conductance of G ~ lGo> as for Pt-H-H- 
Pt contacts typically display step- like decreases in 
dl/dV when conduction electrons are backscattered due 
to energy loss in exciting a local vibration mode Q ■ Fig- 
ure EJc) shows an example of a regular dl/dV spectrum 
for a Pt-D-D-Pt contact. The I-V curve of Fig. [TJb) 
shows clear fluctuations in the voltage window in which 
the peaks in dl/dV appear. At high voltages the fluc- 
tuations have disappeared and the junction displays a 
slightly higher conductance. 

We have evidence that the anomalies are related to 
vibration modes of the single molecule. Firstly, for the 
Pt-Hs system the spikes in dl/dV typically appear at the 
same energies as the transverse vibration mode energies 
of a hydrogen molecule H . Figures 0a) to (c) show a 
sequence of spectra taken when a Pt-D-D-Pt contact was 
slightly stretched over 0.2 A. The small bumps on the vi- 
bration mode shoulders in Fig. Ufa) at ± 48 mV evolve 
into full peaks in Fig. |2Jb) and subsequently into a regu- 
lar vibration mode spectrum for a deuterium molecule 
in between Pt contacts in Fig. Etc). Additionally we 
have introduced hydrogen to atomic gold contacts and 
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FIG. 1: (Color online) (a) dl/dV spectrum on a CO molecule 
bridging a Pt contact displaying symmetric positive spikes. 
The red curve shows a fit by the model discussed in the text 
0; (b) I-V measurement on a CO molecule bridging a Pt 
contact 
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FIG. 2: (a) to (c) display the gradual change of dl/dV spectra 
of a Pt-D-D-Pt contact; (b) and (c) show spectra taken for 
the same junction as (a) after stretching by 0.1 A and 0.2 
A, respectively. Histogram (d) of peak positions in spectra of 
junctions for Au-Eb (hatched bars) and AU-D2 (filled bars). 
All data were obtained at T = 4.2K. 



subsequently measured dl/dV. The spectra show simi- 
lar anomalous spikes as for Pt-H2 . The positions of the 
spikes observed in dl/dV spectra for AU-H2 and A11-D2 
are shown as histograms in Fig. d) • It can be clearly 
seen that the histogram maximum for Au-D 2 at 29meV, 
has shifted by the square root of the mass ratio from 
42meV for Au-H 2 . From the evidence in Fig. [21 it can 
be concluded that the observed singularities are associ- 
ated with the molecular vibration mode energies in the 
contact. 

Spikes in dl/dV point to abrupt changes in the con- 
ductance of the contacts in small voltage windows. Since 
similar spikes have been observed for different molecules 
this suggests that it involves a general feature of single 
molecule junctions. The 'noise' visible in Fig. QJb) is a 
result of two-level fluctuations, that are excited at high 
bias. However, they are apparently not excited directly 
in view of the connection with vibration mode energies 
shown in Fig. [21 Reference Q has shown that by exciting 
internal vibrations of a pyrrolidine molecule on a metal 
surface, the molecule can flip to a different conformation. 
The different conformations have different conductances 
leading to a dip in dl/dV. This forms a basis for our 
proposed explanation. 

We attribute the abrupt changes in contact conduc- 
tance to abrupt switching between two slightly different 
local geometrical configurations that have different con- 
ductances. Current induced two-level systems in meso- 
scopic contacts have been studied in detail in the past 
(see e.g. U 03 EH). In those cases spikes are only 
observed in the second derivative d 2 I/dV 2 , in contrast 
to spikes in dl/dV in our measurements. Since the 
anomalous spikes appear around the vibration mode en- 
ergies, we assume that the switching is induced by intra- 
molecular vibrations. STM studies on single molecules 
on metal surfaces have clearly shown, e.g., that by excit- 



ing an internal vibration mode of a molecule the lateral 
hopping rate of a CO molecule 01 1 * ne rotation of a O2 
molecule or of a C 2 H 2 molecule can be increased 
by two orders of magnitude. In our case the molecules are 
chemically linked by two leads preventing the molecules 
of jumping away and the current is much higher, leading 
to orders of magnitude higher rates. 

We derive a simple model of a vibrationally induced 
two-level system (VITLS) that reproduces the shape of 
the anomalous spikes in our dl/dV spectra very well. 
The model is schematically depicted in Fig. [31 The 
molecule experiences a double well potential landscape 
with ground states 'Frj and <Fi in the two energy minima. 
The two energy minima are separated by a large potential 
barrier that prevents tunneling between the two ground 
states. The molecule in the contact can be vibrationally 
excited by conduction electrons to the state ^> e . In order 
to displace the molecule, the barrier for the reaction coor- 
dinate has to be overcome. The direction of the internal 
vibration needs to be in the direction of the reaction co- 
ordinate (i.e the movement) of the molecule [l^. In the 
case of the Pt-D 2 junctions the spikes have only been ob- 
served at energies corresponding to vibrations perpendic- 
ular to the current direction (i.e. transverse modes). The 
longitudinal vibration mode at about 90 meV [j| never 
showed anomalous spikes. In our quasi one-dimensional 
junctions the transverse modes are the ones that have 
the proper reaction coordinates to induce changes in the 
geometry of the junctions. 

The expression for the current passing through the 
junction can be approximately written as: 
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FIG. 3: Energy landscape of a molecule in an atomic contact. 
The molecule can be vibrationally excited and relax into the 
other energy minimum. 
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We assume that the conductance cr of well and a\ of 
well 1 are different. A step in I(V) and thus a peak in 
dl/dV can be obtained once *i becomes populated. The 
total current is determined by the sum of the products 
of the occupation numbers with their respective conduc- 
tances. For simplicity we take the average of ao and <j\ 
to be the conductance of the excited state. In order to 
determine the shape of the dl/dV spectrum, the depen- 
dence of no, n\ and n e on V has to be determined. Let us 
therefore consider dynamic equilibrium of the occupation 
numbers of all levels involved in the system as displayed 
in Fig. [3] in terms of the transition rates, 

^ = -w e0 + W Qe = 
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dt 
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Here Wfi represent the transition rates from states i to 
states /. The energy difference between the two ground 
states (Fig. |3J| we take to be Ao- The vibration mode 
energy hu> v is much larger than Ao and kgT < Aq to keep 
the system well defined in the ground state for low bias 
voltages. The transition rates W can now be determined. 
Consider as an example the rate W e Q\ 

W e0 = n J 1 f{E,eV)[\- f{E-huj v ,eV)]dE (3) 

— oo 

We assumed here that the electron density of states is flat 
over the energy window under consideration. Where 7 is 
the electron-phonon coupling and f(E, eV) are the non- 
equilibrium electron distribution functions. All other 
rates can be expressed in a similar way. After consider- 
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ing the low temperature limit, 

limiting to positive biases only we can solve the integrals 
of Eq.El 

J f(E,eV) [1 -f(E + AE, eV)] dE 
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With the results of Eq. we can calculate the rates of 
Eq. J2I . The set of linear equations can be solved an- 
alytically and the resulting occupation numbers can be 
inserted in the expression for the current Q). FigureUJa) 
shows a measurement and fit for a H2 molecule in a Pt 
contact, while Fig.0{b) shows the evolution of the occu- 
pation numbers, which clearly shows a sharp drop in no 
and rise in ni once the vibration mode energy is reached. 
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FIG. 4: (Color online) Several fits of VITLS for single 
molecule junction contacts and atomic chains: (a) H2 in a 
Pt contact; (b) Evolution of occupation numbers no (red), ni 
(blue), n e (magenta) for the fit of (a); (c) H2 in a Au atomic 
chain; (d) H2 in Au atomic chain; (e) O2 in a Ni contact. 
The experimental temperature was 5 K in all cases except in 
(d) where the contact was heated to 20 K. For fit parameters 
sec 6|. 



Figure 0Jc) and (d) show spectra and fits for H 2 in a 
Au atomic chain at different temperatures. In Fig. E{d) 
the regular phonon mode of the gold atoms in the chain 
at ±15meV can be clearly resolved next to the VITLS 
spikes. Figure^c) shows clear negative spikes even when 
large conductance fluctuations are present. These con- 
ductance fluctuations arise from interference of electron 
waves that scatter on defects and impurities close to the 
contacts . Often, the regular vibration mode signal is 
masked by these fluctuations for conductances < I Go, 
but the spikes clearly show. Figure^Je) shows a spectrum 
and fit for O2 in a Ni contact, and a fit was also added 
in Fig. n^a), indicating that the model is not limited to 
hydrogen. The energy splitting between the two ground 
states ^0 an d ^1 is l-3meV and is clearly larger than the 
fitted temperatures, which are about 0.3-0.6meV (~ 4- 
7K) in all cases except for the spectrum of Fig.^fd), for 
which the fitted temperature is I.6meV (^20K) in agree- 
ment with the bath temperature. This implies that for 
V<CV m the system is almost entirely in state ^0 (see 
Fig. HTbW which is a prerequisite for obtaining a suffi- 
ciently large decrease of no and increase of ni when V 
reaches V m . 

The model described above involves vibrational excita- 
tion above the barrier. In the case of hydrogen it is also 
possible that the molecule tunnels through the reduced 
effective barrier after excitation. We have elaborated a 
similar model for the case involving a weak tunnel cou- 
pling between excited states in the potential wells, which 
leads to the emergence of two superpositions of the ex- 
cited states [l7(. The resulting fits are not qualitatively 
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different. Furthermore, we have investigated the possi- 
bility of double excitations in the potential wells within 
the framework of our model. The resulting fits could not 
reproduce the experimental data, due to a much slower 
population of the states of well 1. 

As illustrated by Fig. [21 the conditions for obtaining 
a measurable VITLS signal can be obtained by control 
of the MCBJ: It is possible by pulling and pushing the 
contact very gently to explore the potential landscape for 
the molecule so that the proper conditions for a VITLS 
can be obtained. 

Recently STM measurements on a hydrogen covered 
Cu(lll) surface yielded very similar anomalous dl/dV 
spikes 18] . The mechanism is probably related but 
not identical to ours because the peak positions showed 
strong dependence on hydrogen coverage of the surface. 
Moreover, the observed two-level fluctuations were ob- 
served to remain in the o\ -state above the threshold volt- 
age. Possibly the geometry of the experiment allows for 
the molecules to be excited into positions away from the 
STM tip, after which they are slow to return. 

The IETS measurements by Wang et al. on a 
self- assembled monolayer of octanedithiols have been dis- 
cussed jliJ] in view of the unexpected shape of the fea- 
tures in d 2 I/dV 2 . Several of the stronger features have 
the shape of peaks in the first derivative, dl/dV, but still 
at the energies of known vibration mode energies for the 
molecules. Our model produces a good fit to the data 
by Wang et al. but the temperature in the fit is much 
higher than the experimental temperature. This may be 
attributed to inhomogeneous broadening in the ensemble 
of molecules measured. 

We conclude that the observed spikes in dl/dV mea- 
surements on single molecule junctions are most proba- 
bly due to vibrationally induced two-level systems. Our 
simple model fits the data very well. It will be interest- 
ing to investigate potential energy landscapes for single 
molecule junctions in atomistic models in order to ob- 
tain a more concrete picture of the states involved. The 
VITLS acts as an intrinsic amplification mechanism pro- 
viding a strong spectroscopic signature for local vibration 
modes and the anomalous spikes in dl/dV are expected 
to be a general feature of single molecule junctions. It 
can therefore be exploited as a new spectroscopic tool in 
analog y to the "action spectrum" proposed by Komeda 
et. al. |L2j. Specifically, Au atomic chains with H2, that 
very rarely show regular vibration mode spectra, have a 
vibration mode around 42 meV, as seen from Fig.[2Id). 
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